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The interaction between charged bodies in an ionic solution is a general problem in colloid physics and 
becomes a central topic in the study of biological systems where the electrostatic interaction between proteins, 
nucleic acids, membranes is involved. This problem is often described starting from the simple one-dimensional 
model of two parallel charged plates. Several different approaches to this problem exist, focusing on different 
features. In many cases, an intuitive expression of the pressure exerted on the plates is proposed, which includes 
an electrostatic plus an osmotic contribution. We present an explicit and self-consistent derivation of this formula 
for the general case of any charge densities on the plates and any salt solution, obtained in the framework of the 
Poisson-Boltzmann theory. We also show that, depending on external constraints, the correct thermodynamic 
potential can differ from the usual PB free energy. The resulting expression predicts, for asymmetric, oppositely 
charged plates, the existence of a non trivial equilibrium position with the plates separated by a finite distance. 
It is therefore crucial, in order to study the kinetic stability of the corresponding energy minimum, to obtain its 
explicit dependence on the plates charge densities and on the ion concentration. An analytic expression for the 
position and value of the corresponding energy minimum has been derived in 1975 by Ohshima [Ohshima H., 
Colloid and Polymer Sci. 253, 150-157 (1975)] but, surprisingly, this important result seems to be overlooked 
today. We retrieve the expressions obtained by Ohshima in a simpler formalism, more familiar to the physics 
community, and give a physical interpretation of the observed behavior. 



I. INTRODUCTION 

Poisson-Boltzmann theory is a statistical mean field theory 
that characterizes coarse-grained quantities such as the aver- 
age particle distribution function and the electrostatic poten- 
tial together with thermodynamic variables in systems com- 
posed of many charged and point like particles at thermal 
equilibrium. Despite the technical advances in the dilute 
and strong coupling regime Ol 0, Htl, the statistical model- 
ing of real solutions - often in an intermediate regime - is 
still an open problem |@[]. The PB approximation remains 
a good reference theory for describing the essential features 
of electrolyte solutions at thermal equilibrium. It allows to 
model plasmas in the equilibrium regime, colloidal suspen- 
sions through the famous Cell Model [5], or poly electrolytes 
in solution. Moreover, the increasing interest for the biologi- 
cal mechanisms at the sub-cellular scale leads the community 
to deal with the electrostatic interaction of biological objects 
in solution, as for the case of protein-protein interaction |a], 
protein-DNA interaction [7], DNA-membrane interaction fl, 
etc. 

In the case where one is interested in the effective inter- 
action between two charged bodies surrounded by mobile 
charges, it is frequently useful, given the difficulty of the equa- 
tions that have to be solved, to rely on a one dimensional 
problem to capture the physics of the system |9]. This es- 
sentially amounts to focus on the interaction between two 
parallel charged plates in solution. Besides, approximated 
methods have been developed in the past century to correct 
the ID problem as to take into account the geometric ef- 
fects in the interaction of two mesoscopic bodies, thus in- 
crea sing all the more the interest of one dimensional models 

ililliilliililH]. 

In general, the main quantities to be derived in the one di- 
mensional case are ( i) an expression for the free energy of the 



system in the framework of the Poisson-Boltzmann approxi- 
mation, ( it) a differential equation for the mean electrostatic 
potential and, in order to evaluate the actual interaction be- 
tween the two plates, ( Hi) an explicit expression for the pres- 
sure exerted on each surface. 

Various derivations of the Poisson-Boltzmann approxima- 
tion actually exist. A good review of many ways to obtain 
the Poisson-Boltzmann equation has been presented by Lau 
[Ts] including a saddle point approximation in a path integral 
formulation (see also |4]). Less straightforward derivations 
are also available via the Density Functional Theory (DFT) 
[H^ or exact equations hierarchy [ITJ. Finally a less formal 
procedure has been proposed by Deserno et al, in the field 
of colloid physics, to obtain mean field quantities for charged 
systems [5J. 

Most presentations, despite their different approaches, lead 
to a same formula for the pressure, which amounts to the sum 
of a purely electrostatic plus a purely osmotic contribution. 
One merit of the Poisson-Boltzmann approximation is indeed 
that this formula exactly matches the boundary-density the- 
orem at the Wigner-Seitz cell boundary L 1 8,1 as well as the 
contact value theorem on the charged plates ifTgll . The first 
question addressed in this paper is thus whether or not this in- 
tuitive expression for the inter-plate pressure can be directly 
and exactly derived from the Poisson-Boltzmann free energy, 
without need for additional arguments and for any boundary 
conditions. After having introduced the system and its Pois- 
son Boltzmann free energy in Section [III we derive in Sec- 
tion|III]the expected expression for the pressure and show that 
a particular caution should be taken in the choice of the right 
statistical ensemble when different "external" constraints are 
imposed to the plates, as e.g. at constant potential or at con- 
stant charge conditions. 

The pressure formula predicts the presence of a non triv- 
ial equilibrium distance for plates of opposite and asymmet- 
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ric charge densities. This has been shown in the pioneering 
work of Parsegian and Gingell [9i] who used the linear Debye- 
Hiikel theory in the case of high salt concentrations, and more 
recently, by Lau and Pincus |20] in the framework of the non- 
linear Poisson-Boltzmann equation restricted to the case of no 
added salt. 

The consequences of such an equilibrium on the effective 
behaviour of charged bodies in solution can only be assessed 
by a study of the corresponding energy profile, i.e. a com- 
parison of the energy well depth to the thermal energy. If 
the energy gain at the minimum is small with respect to , 
the two charged bodies will not stabilize in the bound com- 
plex and will behave as in the absence of electrostatic inter- 
action. Quite surprisingly, this aspect of the problem is rarely 
addressed in the contemporary literature. Some authors [211] 
discuss in details how the equilibrium distance (the limit be- 
tween attraction and repulsion) depends on the plate charges 
and on the salt conditions, but do not address the question 
of the depth of the free energy well. Nevertheless, very nice 
analytic expressions for both the position and the energy val- 
ues at the ecjuilibrium position have been obtained in 1975 by 
Ohshima [22!]. The paper by Ohshima deals with the more 
complex case of two parallel plates of given thickness and di- 
electric constant, thus leading to a rather complex notation. 
Nonetheless, the important results of Ref. i22ll are worth be- 
ing reproduced today at least in the more usual case of two 
charged surfaces, in that they represent an exact and synthetic 
description of their interaction whatever their charges and the 
ionic strength of the solution. 

In order to illustrate the system behavior in the simple but 
crucial case of monovalent solutions, in Section |IV] we first 
solve explicitly the Poisson-Boltzmann problem and obtain 
the pressure and energy profiles. Then, we focus on the origin 
of the energy minimum and derive an expression for its po- 
sition and depth in the framework of the Poisson-Boltzmann 
theory. We check the agreement between the analytic expres- 
sion and the behaviour obtained by direct numerical integra- 
tion of the Poisson-Boltzmann equation. Finally, we discuss 
the physical origin of the results by investigating the role of 
the different parameters, as the plate charges and the salt ions 
and counter- ions. 



n. THE POISSON-BOLTZMANN FREE ENERGY OF THE 
TWO PLATES SYSTEM 

We are interested in the thermodynamic properties of a sys- 
tem composed of a fixed distribution of charges and of 
point-like mobile ions in a solution at temperature T . The va- 
lence, mass, position and momentum of the ion indexed by "/" 
are denoted by z,, m,-, r,- and /?,, respectively. The Hamiltonian 
of the system can be written as follows: 



^^(W,{p}) = Hk,„+H^ 



pot 



(1) 



where a is the fixed volumic charge distribution in unit of the 
elementary charge e and £ = £08^ is the dielectric constant of 



the solvent. The function (|)(r) is the electrostatic potential, 
^ - o(r')e 



' ^^4%e\r-rj\ J 47te|r-r'| 



a=l 



47i£|r — r' 



where we introduced the ion density of the species a, na(r) = 




ION RESERVOIR 

FIG. 1: A schematic view of the system considered throughout the 
paper. The two semi-infinite planes of charge density Oq and O/,, 
positioned respectively at j = and x = L, are immersed in the ionic 
solution. An ion reservoir freely exchanges ions with the system. 



We will then consider a system composed of two uniformly 
charged plates separated by a distance L and the electrolyte 
solution between them (see Fig.[T]). The fixed charge distribu- 
tion is then 



o(x) = Oo8(x) + Ol8(x - L) 



(3) 



where Oo and Ol correspond to the surface charge densities of 
the two plates positioned respectively at x = and x = L. The 
system is in contact with an infinite salt reservoir Only the x 
coordinate is relevant due to the translation invariance along 
the z and y directions. In the following, we will focus on the 
volume delimited by a given finite surface A of the two facing 
plates. 

As usual, we can obtain the free energy of the system from 
the system partition function Z, as F = —ksTlnZ. The ki- 
netic part Zkin can be easily calculated ll23ll and reads Zki„ = 



U'Lii^a^^'VNal), where Aa = h (IwiiakBT)- ' ' ' is the de 
Broglie thermal wavelength. The potential part of the parti- 
tion function is not that simple to compute, because the elec- 
trostatic part of the Hamiltonian is a function of the position of 
all ions and cannot reduce to a product of uncorrelated func- 
tions. The simplest method to solve the problem is to rely 
on a mean field approximation. The Gibbs-Bogoliubov in- 
equality allows one to find an upper bound for the Helmoltz 
free energy from an average of the Hamiltonian with a trial 
distribution Po(x) plus a Shannon type entropy built from the 
same distribution Po(x). For a given surface A, one gets there- 
fore the following expression for the free energy functional 
per unit surface - that is the Poisson Boltzmann free energy 



x-l/2 
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functional: 



= lim i / -0 

A h^o 2 J- 



(4) 



+ kBlf^ r^'{«OW(ln(A^„0(x))-l)}^/x, 



where we have introduced n^{x) = NaPo{x) following the nor- 
malization relation Na=Aj n1^{x)dx, and the global charge 
density p" defined by 



a=l 



We recognize, in the first term of this functional, the elec- 
trostatic part of the energy of the system, while the second 
term corresponds to the entropic contribution of an ideal gas 
of ions. 

We should therefore minimize the functional Fps [{«"}] /A 
with respect to the relevant functions n^. In order to take into 
account properly the boundaries at x = and x = L we in- 
troduced a parameter h for the calculations and then take the 
limit for h^O. The minimization should be performed under 
the condition of conservation of the whole number of ions of 
type a in the system: this leads to define a generalized free 
energy functional per unit area. 



PB 



f>fi [{«"}] 



P rL+h 
lim TfJa I 



\{x)dx (6) 



where /Ja is a Lagrange multiplier that corresponds to the elec- 
trochemical potential [29] of the ion type a in the reservoir 

The minimization leads to the following relation between 
the mean field ion distributions minimizing fps and the cor- 
responding mean field potential (|)(x): 



(7) 



The reader will recognize in this result an explicit expres- 
sion of the Boltzmann law, here rigorously re-obtained in the 
framework of the mean field approach. 

Together with Eq. (|2|i, giving the electric field as a function 
of the charge distribution in the system, the previous Equa- 
tion (|7]i constitute the solution of the problem. Eq. (|7]i allows 
to obtain a simpler expression for the potential in terms 
of the free and fixed charge distributions in the system. Re- 
calling that the electric potential and the charge density are 
linked by the Poisson equation, i.e. 



Ac^«(r) = -P!i^ 



(8) 



and combining with Eq. we obtain indeed an ordinary dif- 
ferential equation for the adimensional mean field potential 
vj7(ji;) — Pe(|)(x). The resulting fo/iion-BoZfzOTfln« (PB) equa- 
tion reads in our one-dimensional case: 

d^Ml(x) 

\ = -471:4 £ Zana{x), x G \im\+h,L-h] (9) 



dx^ 



with the boundary conditions 



lim^ 

h^O dx 

lim^ 

h^Q dx 



+h 



L-h 



(10) 



where £b = e^/4-nekBT denotes the Bjerrum length and where 
we used electroneutrality of the considered system toget the 
boundary condition. 

Before going on, we should note that it is possible to get 
(5) a constant of motion C by multiplying (|9]l by ^ and then 



integrating in the [+h,L — h] range. One gets 



(11) 



This result will have a crucial role in the definition of the pres- 
sure between the two plates as we will see in the next section. 

In order to compute a general expression for the pressure 
from a thermodynamic definition, we have now to evaluate 
the PB functional free energy at n'^{x) = na{x). The result, 
written in an equivalent but more practical form, will be iden- 
tified with the system free energy. Indeed, at the h ^ Q limit, 
the integrals of any non diverging function in the two exter- 
nal regions vanish. Thus, after an integration by parts for the 
electrostatic contribution, we obtain for the PB free energy 
expressed in terms of the adimensional field \(7: 



P/™[W,a] = hm{J^£;" 



'^Ydx 
dx J 



(12) 



p 

E 

a=l 



L-h 



na{x) [in (A^na) - 1 - pA'a] dx} 



where we used Eqs. ( fTOl i. 



III. DETERMINATION OF THE PRESSURE 

We can now address the problem of finding an explicit 
expression for the pressure 11 on the plates, according to 
the Poisson-Boltzmann theory. In the case of two constant 
charged densities on the plates, i.e. two plates whose charge 
densities are fixed once for ever, the plates self energy is inde- 
pendent of L and the usual derivation of the pressure from the 
free energy can therefore be used: 



n= -A 



dfpB[{na}] 



dV 



3/pB[{na}] 



dL 



(13) 



To facilitate the calculation, let introduce the adimensional 
variable ^ ~ x/L, ifTsIl and define a rescaled electric field E 
for the adimensional potential, £ = — ^ . We then get directly 
from © 



P _ I 

a=l 
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4niB L dt, 



(14) 
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The entropic part of ( fT2] l can be written as 

-1-/1 / TT7 1 ;)<r P 



spB[{n},a] = lim 



by using Eq. (|7]i. We can now take the derivative of Eq. (fTsT l 
with respect to L, and get 



8L 



s/.b[{«},o] = lim / 



1-/! 



1 



d\\td'E 



d d'E 



4714 \dL d^, ^dL 3^ 



a=l 



(16) 



Finally, using again Eq. (|7]i and recalling that ria are func- 
tions of X = t,L, the derivative of the entropic term in the free 
energy, Eq. ( fTST l. becomes 



rl-h 

lim / 



In the same way, we calculate the partial derivative of the 
electrostatic part of Eq. ( IT2b . Grouping the previous results 
and using Eq. ( fTTT i. one gets for the pressure: 



n = C+ hm 

B h^o AkIb 



4k£ 



(17) 



In the present case of constant charge densities on the plates, 
the electric field at the boundaries is independent of L and then 
the second term of Eq. ( fTTb vanishes. The final expression for 
the pressure is therefore: 

n=-— — C=-— — — +kBTl^na{x) (18) 
471:4 87i:^B V "^-^^ / 

The latter result is quite intuitive since it represents the sum 
of the electrostatic stress and the osmotic pressure. It is in- 
deed widely used in the literature. However, in the general 
case of non constant plate charge densities, the second term 
in Eq. ( [TtT i is a priori non zero. Such a term arises for in- 
stance when the potential on the plates is kept constant. Feyn- 
man [24] already pointed out this issue in it's famous course 
on electromagnetism for the case of the pressure between the 
two parallel plates of a capacitor. If, for a given distance L 
between the plates and a given potential difference, we try to 
evaluate the pressure by differentiating the energy of the ca- 
pacitor CU^ /2 with respect to L, then we get two different 
results whether the differentiation is done at constant charge 
or at constant potential. Of course the two derivations should 
give the same result since they refer to the same state of the 
system. The explanation of this apparent paradox is actually 
quite simple: when differentiating at fixed potential we in- 
clude implicitly the energy supplied by the generator to keep 
the potential difference constant while varying L. Since this 
work only modifies the self energy of the plates and not their 
interaction, we have to subtract this part from the result. We 



retrieve in this way the correct result of the constant charge 
case. 

In our case, one can easily realize that, at the boundaries, 
the expression \j73£ /dL is equivalent to \(73a/3L and thus 
corresponds to the energy used to bring charges to the plates 
when varying the distance L. Then, the second term that arises 
in Eq.(fT7]l is exactly the analogous of the generator term in 
the capacitor problem, and appears only because, using Gibbs 
terminology, we don't work in the correct statistical ensem- 
ble. Indeed Eq.(fT3]l represents a very useful and systematic 
procedure to get an expression for the pressure provided that 
the ensemble or equivalently the effective potential is chosen 
carefully. For instance, in the case of constant potential on the 
plates we have to make a functional Legendre transform to get 
the relevant thermodynamic potential: 

_ rL+h _ 

m[{n},Mfe] =/ra[{n},o] -lim / y^{x)o{x)dx (19) 

where {^}b denotes the electrostatic potential at the bound- 
aries. In this case Eq. ( flj] ) must be replaced by: 



dL 



{¥}(, 



(20) 



Since the derivative of the additional term in gpg exactly bal- 
ances the second term of Eq. ( fT7] i. we finally retrieve the gen- 
eral resul t dTSl) . For a different but equivalent discussion, see 
also Ref.iZSll. 



IV. EXISTENCE AND CHARACTERIZATION OF THE 
ENERGY MINIMUM FOR ASYMMETRIC CHARGED 
PLATES 

A. Zero pressure distance Lm,„ 

It can be interesting to illustrate the implications of Eq. ( fTSl ) 
for the case of two charged bodies in a 1:1 solution. The ex- 
tension to a multivalent solution is straightforward |30], but 
the simpler case is more instructive. First of all, we note that 
the expression Eq. ( fTSl l represents only the pressure due to the 
inhomogeneous electrolyte solution between the plates and 
we have therefore to add the pressure contribution that comes 
from the homogeneous electrolyte solution surrounding the 
system. Note that, since it has been assumed that our system 
is electrically neutral, this osmotic contribution is homoge- 
neous in the surrounding space. Let i = nf,._i = nb be the 
bulk concentrations for the positive and negative monovalent 
ions. We thus introduce the excess pressure P that satisfies 



n(x) - Hoc = 

ksT fd^ 



-kBTY^{na{x) 



-nb) 



(21) 



and vanishes when L increases toward infinity. The reservoir 
ions are here assumed to behave like an ideal gas, coherently 
with the mean field approximation. 
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Now, f is a function of L, that can be obtained by solving 
the PB Equation (|9]l for each L e ]0; +°°[. 

We then numerically analysed the sign of P as a function of 
L and of the plate charge ratios r = Ol/oo. For r > (charges 
of same sign) the plates always repel each other, which is a 
general consequence of the PB theory [9]. For the particular 
case r = — 1 the interaction is instead always attractive. In- 
terestingly, in the more general case of r < 0, and r ^ — I, 
there always exists one and only one equilibrium distance 
Lmin ('") between the plates for which we observe a transition 
between attraction and repulsion (i.e. a vanishing P). The 
transition occurs at a distance that depends on the charge den- 
sities of the plates, and a pronounced repulsion always appears 
at short distances, despite the fact that the plates are oppositely 
charged. 

Such transitions were already predicted in linearized treat- 
ments of the problem in 1972. More recently, the non 
linear case has been reconsidered f2T], although an exact 
derivation of the transition distance as a function of the plates 
charges and salt concentration had akeady been obtained by 
Ohshima |22] in 1975. 

Following the main lines of Ref. ||22|] . it is possible to obtain 
an analytic expression for the equilibrium position explicitly 
dependent on the plate charge densities and on the salt concen- 
tration, for the case of a monovalent solution. We perform this 
calculation explicitly in Appendix A. We obtain the following 
expression for the position of the energy minimum L„„„: 



In 



K\{2 + 



-4) 



K\{2- 



-4) 



(22) 



where we have introduced the Debye length, Xd ~ 
1 1 y/^Kls nb and the adimensional charge densities o'o = 
AtzIbXdOq and o' l = 4-k£bXdOl. A similar expression for 
the distance at which P = is given in Ref. f^ll], Eq. (9) lIsTll . 

In Fig.|2]we compare the previous expression Eq. (l22l l for 
Lmin with the corresponding values directly obtained from the 
numerical solution of the PB equation, for different salt con- 
centrations and charge density ratios. The minimum positions 
Lmin are numerically estimated directly from the energy pro- 
files. As expected, the formula of Eq. ( |22] | exactly agrees with 
the numerical results. 

Two limiting regimes can now be considered. Let introduce 
the Gouy-Chapman lengths for both plates, A-o = 1 1/271^^00 1 
and Xl — \l /2k£b<^l\- In low salt conditions, Xd ^ A-o and 
Xd ^ Xl. As a consequence, the position of the energy mini- 
mum is approximately given by 




FIG. 2: Comparison between the estimation for the position of the 
energy minimum of Eq. ( I22l l (solid lines) and the values of l^min 

ob- 
tained on the basis of the direct resolution of the Poisson-Boltzmann 
model (points). We chose Oo = — 0.05e/nm^. Lengths are given in 
nm, and as functions of the ratio r = Cl/oq, for different salt con- 
centrations: 0.001 M (asterisks), 0.01 M (circles), 0.1 M (squares), 1 
M (diamonds). 



counter-ions of the most charged plate is much larger than the 
salt concentration. In this case, the short range repulsion is 
therefore mainly due to the counter-ions of the most charged 
plate. We stress that the solid and dotted curves in Figure |2] 
also correspond to this low salt regime. 



0.5 
x/l 



FIG. 3: Positive (solid line) and negative (dashed line) ion dis- 
tributions in the inter-plate space, at the equilibrium distance L — 
Lmin, for the case of two charged plates with charge densities Oq = 
— 0.05e/nm^ and ai = 0.1 e/nm^, and for two different salt concen- 
trations nt,: 0.001 M (low salt, a) and 1 M (high salt, b). The thick 
grey line indicates the value of the bulk ion concentration nt,. 



• 2Xd 



1 



1 



lA-L-A-ol (low salt), (23) 



Inversely, at high salt, Xd ^ Xo and Xd <C Xl, and the equi- 
librium is 



In this limit, L„„„ becomes therefore independent of the salt 
concentration and is only a function of the plates charge den- 
sities, i.e. of the ratio r for the cases considered here since A,o 
is always kept fixed. In Fig. [3] we report the local ion distri- 
bution in the inter-plate space as a function of x/L for a given 
choice of the plate charges and for different bulk ion concen- 
trations Hh- At low salt (Fig. [3] a), the concentration of the 



' Xd 



1 1^0 



- Xd 



1 1^0 I 

In — 



(high salt). (24) 



In this limit the estimated equilibrium length Lmin is then pro- 

portional to the Debye length, i.e. to . As shown in Fig. [3] 
b, the short range repulsion is indeed essentially due to the salt 
ions whose osmotic effect is modulated by the charges on the 



6 



plates. The dashed and dot-dashed curves in Figure |2] corre- 
spond to the high salt regime. In this high salt regime, a good 
approximated expression for the equilibrium position can also 
be obtained in the framework of the linearised PB equation, 
as expected. We checked indeed that the resulting expression 
^ matches well the curves in Fig. |2] for any ratios r at high 
salt concentrations (data not shown). Instead, the linear PB 
approximation cannot reproduce the observed behavior at low 
salt, whereas the expression of Eq. (l22T i remains exact. 




0.5 1 1.5 2 2.5 3 3.5 4 



— r=-2 r=-l,2 — r=-l--r=-0.8 — r=-0.5- ■r=0 ■ r=l 



FIG. 4: The interaction energy per unit surface E (ksT/ 100 nvn?) 
for the interaction of a plate of charge density Oo = — 0.05e/nm^ at 
X = with different plates of charge densities Cl, as a function of 
the distance L (nm) between them. The ratio r = Oi/Oo between the 
densities varies from —2 to 1 according to the figure legend. The 
plates are immersed in a 0.1 M monovalent solution. 

How the P = condition should be interpreted in terms of 
electrostatic and osmotic contributions? The mechanism lead- 
ing to an equilibrium position is not difficult to understand, 
starting from the expression of the excess pressure, Eq. ( |2TI ). 
As already observed, the pressure is proportional to the con- 
stant of motion C and is therefore constant in the inter-plate 
space X e [0,L]. Let then consider the pressure P exerted on 
one plate, e.g. at x = 0. In this limit, the electrostatic term 
in Eq. ( 1211 1 simply reads —hgT Oq /2, and is therefore inde- 
pendent of L. On the contrary, the osmotic term depends on 
the ion concentration and is therefore a function of L. We 
then note that the expression for the ion density, Eq. (|7]i can 
be rewritten in terms of the bulk concentrations and of the re- 
duced potential as 

«„(x)=«fce--'«^W. (25) 

The equilibrium condition P = 0, when calculated both on the 
left and right plates, can thus be written as a condition for the 
mean field potential at the boundaries that reads 

cosh%(L,^) = ^^^ + l. (26) 
nb 

cosh^i(L«/) = + 1 . (27) 

nb 

where we introduced Vo(^) — — 0) xj7^(L) = \(7(x = 
L) VL € [0, oo[, and where Leg is a distance between the plates 
for which P = 0. 
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— r=-2 r=-1.2 — r=-l--r=-0.8 — r = -0.5- -1=0 ■ ■ ■ r=l 



FIG. 5: The adimensional mean field potential at x = 0, ^7(0), as 
a function of the inter-plate distance L (nm) and for the same ionic 
strength and charge density ratios of Fig.|4l The two horizontal lines 
correspond to the two roots ± arccosh(7l£gaQ/ii/, + l) of the left 
plate equilibrium condition, Eq. M6\ . The circles emphasize the non 
trivial solutions leading to an equilibrium position. 



From now on we will focus on the case of asymmetri- 
cally and oppositely charged plates, r < 0. In Figure |5] we 
compare the adimensional potential V(7o(L) to the two roots 
±arccosh(7t£BOQ/nf, + l) of Eq. ( |26] |. A trivial solution to 
Eqs. (|26] i and ( |27] i corresponds to the limit L-^ where P — 
by construction. In this limit, and typically for L^Xo + Xl, 
each plate charge is neutralized by its cloud of counterions 
as if the other plate didn't exist. At a large enough distance 
X from the plates, therefore, the ionic atmosphere behaves as 
an ideal gas and its pressure contribution Tl{x) exactly com- 
pensates the reservoir pressure IIoo. We also note that, in this 
case, according to ( |25] ), the sign of \(7q and \(7^ at infinity are 
opposite for oppositely charged plates. 

Let now look for another, finite solution of Eqs. ( 1261 ) and 
(|27| | leading to Leg = L,ni„. For this case, we only have to 
determine the sign of \(7q (L„„„ ) and '^i^{L,nin)- It easy to see 
that the only possible solution is that the counterions of one 
plate prevail on both plates, so that the signs are both equal. 
It is the potential of the less charged plate that will change 
its sign at a given distance L between infinity and L,„,„. In 
Figure |5]we show indeed that the sign of the adimensional 
potential vj/g changes for r < — 1 whereas it remains the same 
when r > — L The opposite arises for ij/^ (data not shown). 

The physical interpretation of the observed equilibrium is 
thus straightforward. The most charged plate carries its cloud 
of condensed counterions when approaching the less charged 
one. The counterions of the less charged plate are more easily 
released in the bulk when the two ion clouds overlap. This 
process continues until the electroneutrality constraint pre- 
cludes a further release of ions, this leading the ion concentra- 
tion, and therefore the repulsive osmotic pressure, to increase. 
The equilibrium is then obtained at the distance L„„„ where 
electrostatic attraction and osmotic repulsion exactly balance. 
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FIG. 6: Comparison between the estimation of the energy minimum 
of Eq. l l28t (solid lines) and the values of the energy at the min- 
imum position L,„,„ obtained by direct integration of the Poisson- 
Boltzmann model, for different salt concentrations: 0.001 M (aster- 
isks), 0.01 M (circles), 0.1 M (squares), 1 M (diamonds). We have 
udes again Oq = — 0.05e/nm^. In order to use more natural units, 
energies are given in units of ksT / 100 nm^. 



imum value for r = — 1 , i.e. when the position of the minimum 
degenerates to L = 0. Energy depths are up to roughly 10 kgT 
per 100 nm^. Figure |6] also confirms that the minimum depth 
becomes constant for r < — 1, and coincides in this case with 
its (maximum) value at the singular value r = — 1 . 

In his paper, Ohshima 1E2I1 also obtains an analytic expres- 
sion for the energy at the minimum. An equivalent calcula- 
tion, adapted to our formalism, is presented in Appendix B. 
The final result reads 



-4 — 2 — |o ml arcsinhl 



(28) 

where &,„ is the adimensional charge parameter related to the 
smallest plate charge density, i.e. oJ„ — min(OQ,0^). In Fig- 
ure |6] we compare the previous expression for the value of the 
energy depth with the results obtained by direct integration of 
the Poisson-Boltzmann equation. The two results show a per- 
fect agreement. Together with Equation (l22T l. the last result 
allows to a rapid and precise estimate of the equilibrium po- 
sition and strength, and represent therefore a powerful tool in 
order to study the effective interaction between charged bod- 
ies in solution. 



B. Energy value £■„„„ at the minimum 

We have shown that the pressure always vanishes at a given 
inter-plate distance for oppositely, asymmetric charged plates. 
Nevertheless, if the presence of a vanishing pressure always 
corresponds in principle to an equilibrium position between 
the plates, it does not guarantee by itself that this equilibrium 
position will be stable enough to be relevant from a thermo- 
dynamic point of view. Indeed, in order to assess the real 
existence of a stable equilibrium, we need to estimate the cor- 
responding energy gain. We stress again that, if the energy 
gain at the minimum is small with respect to ksT, the two 
charged bodies will behave as in the absence of electrostatic 
interaction (the energy going to zero at large distances). On 
the contrary, a deep minimum will make the bodies stabilize 
at a non-zero equilibrium distance. 

The explicit calculation of the whole function P{L) allows 
us to evaluate the energy profile E{L) and compare the depth 
of the potential well to the thermal energy ksT. The ener- 
gies per unit area are shown in Fig. |4] (bottom), for the same 
charge densities as considered above. In order to use more 
natural units, energies are given in units of ksT / 100 nm^. 
As expected, an energy minimum always exists for r < and 
ry^—l. Interestingly, while for < r < 1 the energy minimum 
depth shows a relevant dependence on the second plate charge 
density Ol, this dependence disappears for r < — 1 where the 
depth becomes constant. 

A more systematic investigation of the energy minimum 
depth for varying charge densities and salt conditions is shown 
in Figure |6l The value of the energy per a unit area of 100 
nm^ for different ionic strengths is given as a function of the 
ratio of the plates charge densities r, again for Oo fixed at 
— 0.05e/nm^. In low salt, the energy minimum depends on 
the ratio r and on the salt concentration, and it reaches its max- 




FIG. 7: Sketch for the calculation of the interacting area for the case 
of two identical spheres. The area of each dashed spherical surface 



Nevertheless, we have to stress finally that Equation 
as well as Figures |4]and|6] only give the energy per unit area 
(fixed to 100 nm^ in the Figures). To obtain the total energy 
between two charged bodies in solution and compare it to the 
thermal energy, the surface and geometry of the bodies should 
be taken into account. This roughly amounts to multiply the 
energy by an effective interaction area, but the estimation of 
this area is not easy in that it depends on the bodies shape. In- 
deed, the variation of the interaction with the distance should 
be included in the calculation of the effective interaction area. 
A typic al choice is the use of the Derjaguin approximation 
II271I28II . that calculates the interaction energy U between two 
curved surfaces by integrating the interaction energy per unit 
area between two flat plates E{L) as 



U: 



E{L)dA~f{[Ri],[R2]) E{L)dL, (29) 
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where D,„,„ is the distance of closest approach between the 
two curved surfaces, dA is the differential area of the surfaces 
facing each other, [Ri] and [R2\ represent the sets of principal 
radii of curvature of the surfaces 1 and 2, respectively, at the 
distance of closest approach, and f{[al], [a2] ) is a function of 
the radii of curvature of the surfaces. A very rough estimate 
should consider that the interaction between the two surfaces 
becomes negligible when their distance becomes larger than 
the screening Debye length Xo- In a 0.1 M solution, the De- 
bye length is of the order of 1 nm. For the case of two spheri- 
cal colloids of 100 nm diameter, (it varies from approximately 
0.3 nm for an 1 M solution to 10 nm for 0.01 M). If the two 
spheres are in contact (i.e. for L = 0), a simple geometric con- 
struction (depicted schematically in Fig. IT) allows to calculate 
the surface area on each sphere that is separate by less than a 
Debye length from the facing one, this leading to a surface of 
the order of 300 nm^ (roughly 1 % of the whole sphere sur- 
face). As a consequence, the depth of the energy well would 
reach approximately ^20 kgT in these conditions, and be thus 
large enough to ensure its stability at ordinary thermal condi- 
tions. Note however, that the two spheres will nothe in contact 
anymore at equilibrium, but at a distance close to L„„„, which 
can vary from to approximately 2 nm depending on the value 
of the ratio r. When this distance becomes comparable to A-d, 
the interaction area decreases considerably. Therefore, the ac- 
tual interaction surface will be in general reduced to a value 
that depends on the two lengths Lmm and A-o, and should be 
calculated case by case. Note that, in this sense, the shape of 
the two interacting bodies is bound to play a crucial role, in 
that it can modify considerably the distances between the fac- 
ing surfaces and consequently the effective interaction area. 



derivation potentially leads to underappreciate the role of the 
boundary conditions. 

We then explicitely solved the problem in the simple case of 
a 1 : 1 salt solution and observed a very rich behavior as a func- 
tion of the ratio of the plates charge densities and of the salt 
concentration, with a non trivial equilibrium position arising 
in large intervals of these parameters. The distance at which 
the osmotic and electrostatic pressures are in equilibrium is 
finely tuned by the system parameters. Such equilibrium po- 
sition can stabilize two asymmetrically charged bodies at a 
nonzero distance, provided that the corresponding free energy 
gain is large enough compared to the thermal noise. At a given 
temperature, the stability of the complex can therefore be as- 
sessed only by explicitly calculating the depth of the corre- 
sponding energy well. We obtained a readily available answer 
to this problem by deriving analytic expressions for the posi- 
tion and depth of the energy well, by a reactuahzed version of 
the overlooked derivation of Ohshima (tS] . 

In order to compare the energy values with the thermal 
energy, an estimation of the interaction area is also necessary. 
As an example, we gave an estimation for the typical case 
of spherical colloids, and found that the interaction energy 
well is typically of the order of several ksT, this leading to a 
quite deep minimum. This estimation is only a rough approx- 
imation because, for a given problem either in biological or 
colloidal systems, the behavior of the two interacting bodies 
is strongly dependent not only on salt conditions and on the 
bodies charge but also on their shape. The relevance of the 
minimum is consequently rather delicate to compute and 
model-dependent. 

APPENDICES 



V. SUMMARY AND CONCLUSION 



A. DETERMINATION OF THE MINIMUM POSITION. 



In this paper we derived, within the genuine framework of 
the Poisson-Boltzmann approximation, the main quantities of 
interest for the interaction between two charged surfaces in 
solution. 

Interestingly, such a rigorous derivation brings up an ex- 
tended formula for the interplate pressure, formally differing 
from the usual expression (giving the sum of electrostatic and 
osmotic contributions) by an extra-term. While in the con- 
sidered case of fixed charges boundary conditions, this addi- 
tional term vanishes and we recover the standard expression 
for the pressure, the question arises whether it leads to a mod- 
ified result when different boundary conditions are consid- 
ered. We recalled that the Poisson-Boltzmann free energy is 
properly used as the reference thermodynamic potential only 
at fixed charges. To show how the choice of the thermody- 
namic potential depends on the boundary conditions, we ex- 
plicitely solved the problem in the fixed potential case, and 
show how to recover the physical meaningful pressure in that 
case. These precisions on the interplay between boundary 
conditions and thermodynamic ensemble have been obtained 
here thanks to a detailed and deductive derivation from the 
very basis of the Poisson-Boltzmann approach. Let us stress 
once again that shortcutting the details of the mathematical 



Starting from Eq. ( fT2l l and taking into account the contri- 
bution of the surrounding ions as in Eq. ( |2TI ). one can obtain 
a suitable form for the excess free energy that will allow one 
to evaluate its amplitude at a given position x. We start by 
writing the excess free energy as 



dx 



2 

-2nhL+ ^ 
a=l 



-,70 



na[-ZaV(/(x)] dx (Al) 



where the 2ni,L term arises from the Ha, contribution. Now, 
using (fTTT i and (|9|l one finds: 



pJpB[o{n}] = 2«JQc'+l ]L 



\|/(jc)sinhi|/(^) dx 



2nh ■ 



-I n"^ 



dx 



dx 



(A2) 
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where we have introduced the Debye length, Xd = 
1 / \/87i£b lib, C = Cj {Anisnb) and we performed an integra- 
tion by part for the last equation. 

In order to calculate L,„,„, it is convenient to introduce the 
following adimensional parameters: 



ri = x/Xd 
Y(ri)=\(7(x) 
dy 



9(11) 



(A3) 
(A4) 

(A5) 



In terms of these variables, the Poisson-Boltzmann equation 
just writes 



— — smhY 
dy\ 



(A6) 



with the boundary conditions 



9(0)- 

e(iiL) 



-o'o 



- 471^6 A.£,Oo 



(A7) 



and the equivalent of Eq. ( fTTT l. giving the constant of motion 
for the system, reads 0^ = 2coshY+C'. From this equation 
one gets 



(02 



= cosh^ y = 1 + sinh^ Y 



4 

(02-C')2 



-1 = 



dQ 
dr[ 



dr[ 
dQ 



1 



(A8) 



v/(02-C')V4-l 
By integrating r| from to r[L, with |(ir| | = dr[, we thus obtain 



dd 



We are interested here in the case of oppositely charged plates, 
where the energy minimum does exist. We have therefore 
to consider different cases. Let assume for the moment that 
|o'o| > |o'l|. Ifnowo'o <Oando'L >O,then0(O) and0(riL) 
have both positive values and |ii0| = ~dQ. We then have 



/■9(i1l) 


dr[ 


r 

dQ^ 


/e(o) 


dQ 


Jo 



^L = 





c/ri 




dQ 



dQ 



o'o<0, o't>0. 



(A9) 



On the other hand, if o'o > and o'l < 0, then 0(0) and 0(r|L) 
have both negative values and \dQ\ = dQ. We thus obtain 





dr] 




dQ 



dQ 



o'o>0, o'l<0. (AlO) 



Nevertheless, since | dr\/dQ\ is only a function of 0^ we can 
make the change of variable © = —0 in Eq. dAlOb and retrieve 
the result iA9i . Therefore, we always have 



o'lI ^(02_C' + 2)(02-C'-2) 



dQ (All) 



If now L = L„in, then P = 0, i.e. 11 = Yloo = IksTrib, and 
therefore C' — —2. At the equilibrium position, Eq. dAl II ) be- 
comes then 



0v^ 



dQ 



Ic^'ol/Z 



1 



o'z.l/2 avc?+T 



-.dQ 



where we introduced a = 0/2. As a primitive function of 
1 / {x\Jy?- + 1) is In [xj (1 + \J yp- + 1)^, we get for the case 
lo'ol > Io'lI 



:lnl 



|o'o|(2+^^^rF+4) 
|o'l|(2+v/|o'oP+4) 



(A12) 



The extension to the opposite case of |o'o| < |o'l| is 
straightforward, this leading to the following final formula for 
the position of the energy minimum L,„,„ = A-Dr|i,,„„ : 



In 



Kl(2- 



.'2 



-4) 



Kl(2- 



-4) 



B. DETERMINATION OF THE ENERGY VALUE AT THE 
MINIMUM. 

Following the main lines of the calculation presented in 
ll22il . we here look for an analytic expression for the energy at 
the minimum. In our framework, the interaction energy com- 
puted numerically directly from the integration of the excess 
pressure writes formally (by definition of the integration): 



p£(L) = P(j(L)-iF(oo)) 



(Bl) 



where J (L) is given by ( IA2l l. By using the dimensionless 
parameters defined in the previous section, we have 



P^F = {{K{L)+J(L)-l{L)}lnb'ku, 



(B2) 



with 



J{L) 

m 



(\c' + i)y\L. 

Y(il)0(ri)J^ ., 

■TlL 

Q^dy\. 





We will now, first, calculate the three contributions to 
f {Lmin), and then the corresponding contributions to f (oo). 

Let start by calculating I{L„,i„). From the definition, we 
have 



I{L„ 



liri 
dQ 



20 
02T4 



dQ 



2(V|o'or 



-4) 



dQ 



(B3) 
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where we used the fact that, on the [0,r|L] range, d^/dQ = 
-\d^/dQ\. 

To calculate J{L„,i„) we recall that when L — Lmi,, we have 
coshY(O) =a'o/2 + l. By using the relation (coshx— l)/2 = 
sinh^(jc/2), we then get sinh^ y(0) /2 = o'o/4, and thus 



(B4) 



(B5) 



|sinh^|=.i|o'o| 
In the same way we have 

|sinhI^^|-i|o'i|. 
We can now easily calculate 7(L„„„) starting from 

JiL„,i„) = Y(iiL,„„,)e(iiL,„„,)-Y(o)e(o) 

= Y(ilL,„„,)|o'L|+Y(0)a'o. 



Now, at L,„i„, 'Y{^L„,i„) and y(0) have the same sign, which is 
governed by the most charged plate. Let again focus on the 
case where o'o < and |o'o| > \o'l\, as used in our illustra- 
tions. In this case, the sign of y(i1l,„,„) and y(0) is the same as 
the sign of a'o (cf. Figure 4), this leading to 

7(L,„,„) - -2|o'i| arcsinhl^l -2|o'o| arcsinh|^| (B6) 
We easily obtain that K{L,„in) = since C' = —2 for L = 



The overall result for (L„„„ ) reads therefore 
{L,„in) = Anblo S arcsinhl^l 



(B7) 



lo'ol arcsinhl^l -(^|o'o|2+4- ^|o'l|2 + 4; 



Let now calculate J {°°). We have to be quite cautious to 
compute I{°°). Indeed, we have 

and the point is that dT\/dQ has not the same sign all over the 
range [0, Actually because of the infinite distance between 
the two plates, each plate tends to behave as a single plate in 
this limit, and thus there exists a distance for which 9 = 0. 



So, 9 will be initially equal to |a'o|, then decrease to zero and 
increase again to |o'l|. We have therefore: 



9^ 



dT\ 


r\<^'L\ 
dd+ 9^ 


c/rj 


dQ 


Jo 


dQ 



dQ 



2v^ 



2^9^+4 



/(oo) 



= -8+2.^|o'o|2+4 + 2.^|o'l|2+4. 

Besides, we have fory(oo): 

JH = YH|o'L|-Y(0)|a'o| 

= 2|o'l| arcsinhj^^l — 2 |o'o| arcsinh|^^| 



(B8) 



(B9) 



because it's only the reduced potential Y corresponding to the 
plate with the lowest charge (in absolute value) that changes 
its sign between L„„„ and Loo (cf. again Figure 4). 

When there is an infinite distance between the plates we 
also have C' = —2 (intuitively because there is no more inter- 
action between the plates) and thus K{°o) ~ 0. We then have: 



(oo) ^AribXo \ |o'L|arcsinh 



(BIO) 



-|o'o|arcsinh|^|-(-4+^|o'o|2+4- 



o'lP+4) 



From the evaluation of Equation ( IBll ) at L = L„„>,, we get 
then finally the following expression for the energy at the min- 
imum in the case when o'o < and lo'ol > |cj'l|: 



p£'mm — Sn/jA-D 



Io'lP- 



- 4 — 2 — I o'i I arcsinh 



(BID 

On the other hand, one can easily be convinced that inverting 
the roles of |o'o| and |o'l| leads to the same expression as 
Eq. (IB 11 1) where |o'o| and |o'l| are inverted. Therefore, the 
very general result writes (Eq. (|28] |) 



: 8n/,A,D 



t' 12 _ 



-4-2- 



arcsinhl 



(B12) 



where o'„, is related to the smallest plate charge density, i.e. 
|o;„|=min(|o;,Uoi|). 
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